The switching dynamics of the bacterial flagellar motor 
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Many swimming bacteria are propelled by flagellar motors that stochastically switch between 
the clockwise and counterclockwise rotation direction. While the switching dynamics are one of the 
most important characteristics of flagellar motors, the mechanisms that control switching are poorly 
understood. We present a statistical-mechanical model of the flagellar rotary motor, which consist 
of a number of stator proteins that drive the rotation of a ring of rotor proteins, which in turn 
drives the rotation of a flagellar filament. At the heart of our model is the assumption that the rotor 
protein complex can exist in two conformational states corresponding to the two respective rotation 
directions, and that switching between these states depends on interactions with the stator proteins. 
This naturally couples the switching dynamics to the rotation dynamics, making the switch sensitive 
to torque and speed. Another key element of our model is that after a switching event, it takes 
time for the load to build up, due to polymorphic transitions of the filament. Our model predicts 
that this slow relaxation dynamics of the filament, in combination with the load dependence of the 
switching frequency, leads to a characteristic switching time, in agreement with recent observations. 



I. INTRODUCTION 

The bacterium Escherichia coli can swim toward at- 
tractants and away from various noxious chemicals. It 
is propelled by several flagella. Each flagellum is under 
the action of a rotary motor, which can rotate either in a 
clockwise (CW) or a counterclockwise (CCW) direction. 
When all the motors run in the counterclockwise direc- 
tion, the flagella form a helical bundle and the bacterium 
swims smoothly. When one motor switches direction to 
run in the clockwise direction, however, the connected 
flagellar filament disentangles from the bundle, and the 
bacterium performs a so-called tumble. These tumble 
events randomize the cell's trajectory, and it is the mod- 
ulation of their occurrence that allows E. coli to chcmo- 
tax. Here, we present a statistical-mechanical model that 
describes the switching dynamics of the rotary motor. 

A cartoon of the bacterial flagellar motor is shown in 
Fig. [T] It consists of a protein complex called the 
rotor, and a number of stator proteins that are fixed in 
the inner membrane and the peptidoglycan layer. Inter- 
actions between the stator proteins and a ring of FliG 
proteins of the rotor protein complex drive the rotation 
of the rotor, and thereby the rotation of the flagellum, 
which is connected to the rotor. The rotation direction is 
determined by the concentration of the phosphorylated 
form of the messenger protein CheY, which binds to the 
ring of FliM proteins of the rotor protein complex. 

The bacterial flagellar motor has a number of char- 
acteristics that reflect its design principles. One is the 
torque-speed relationship of the motor, which has re- 
cently been modeled by Oster and coworkers 0]. An- 
other is the power spectrum of the motor switching dy- 
namics, which reflects the time scales on which the motor 
switches direction. Recently, Cluzel and coworkers mea- 
sured power spectra for mutant cells by monitoring the 
rotation of a 0.5 //m latex bead connected to a flagellum. 
S In these cells, a mutant CheY protein, CheY*, 
was stably preexpressed 0, . This protein mimics the 




FIG. 1: Cartoon of the fiagellar motor of the bacterium E. 
coli. The figure is courtesy of DeRosier [J] . 



effect of CheYp, but its concentration is not affected by 
the dynamics of the chemotaxis signaling network. The 
power spectra of these mutant cells therefore reflect the 
intrinsic switching dynamics of the motor. 

Intriguingly, the power spectra of these CheY* mutant 
cells are not consistent with a two-state Poisson process, 
in which the switching events are independent and the 
CW and CCW intervals are uncorrelated and exponen- 
tially distributed They exhibit a distinct peak at 
around Is""'^ 01, which means that there is a character- 
istic frequency at which the motor switches. Cluzel and 
coworkers [1] suggest that an earlier model developed by 
Duke and Bray might be able to explain the power 
spectrum. This model, however, is based on an Ising 
system, which is a mesoscopic equilibrium system. Such 
a system cannot exhibit a peak in the power spectrum 
00,@- The peak means that the switching dynamics is 
coupled to a non-equilibrium process. 

We argue that to explain the switching dynamics of the 
motor, we have to integrate a description of the switch- 
ing dynamics of the rotor with a description of both the 
flagellum dynamics and the dynamics of the stator pro- 
teins that drive the rotation of the rotor. The rotor pro- 
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tein complex is modeled as an Monod-Wyman-Changeux 
(MWC) model 0], in which the proteins of the complex 
collectively switch between a clockwise and a counter- 
clockwise conformational state. Interactions between the 
stator proteins and the FliG proteins of the rotor do not 
only drive the rotation of the rotor, but also continually 
change the relative stability of the two conformational 
states of the rotor. This couples switching to the non- 
equilibrium process of rotation, and makes the switch 
sensitive to torque and speed. Our model predicts that 
the probability for the rotor to switch increases with the 
load in the zero-load regime, as observed experimentally 
(lot . But, to fully describe the switching of the motor, 
the rotor's switching and rotation dynamics have to be 
coupled to the conformational dynamics of the flagellum. 

Bacterial flagella can exist in different polymorphic 
states [nl, m, 13 13 , which are either left-handed or 
right-handed helices. When the motor runs in the CCW 
direction, the flagellum adopts a left-handed, normal 
state [l^, while if it runs in the CW direction, it adopts a 
right-handed, semi-coiled or curly state By pulling 
on a single flagellum using optical tweezers, Darnton and 
Berg recently observed that transitions between these 
polymorphic forms occur in discrete steps, during which 
elastic strain energy is released [13]. We argue that the 
change in the torque upon a motor reversal induces a 
polymorphic transition that proceeds via a similar series 
of discrete steps. Since in each of these steps strain en- 
ergy is released, the torque on the motor, and hence the 
switching propensity, remains low. However, when the 
flagellum ultimately reaches its final polymorphic form, 
the strain energy can no longer be released and the torque 
on the motor increases. Due to the exponential depen- 
dence of the motor's switching propensity on the applied 
torque, the switching propensity now increases rapidly, 
giving rise to the characteristic switching time 
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FIG. 2: Energy surfaces for the interaction between the rotor 
complex and one stator protein. A: Energy surfaces corre- 
sponding to the two conformational states of the stator for a 
given conformational (CW) state of the rotor protein complex, 
according to 0, ■ The dashed lines denote the hopping win- 
dows. The thermodynamic driving force is AG = 2e x pmf. 
B: Two energy surfaces of the motor, corresponding to the 
CCW and CW states of the rotor, for a given conformational 
state of the stator protein; the two surfaces are assumed to 
be each other's mirror image. In total, each rotor-stator in- 
teraction is characterized by 4 surfaces, corresponding to the 
2x2 conformational states of the stator and rotor proteins. 



advanced by an angle of 27r divided by 26, the number of 
FliG proteins within the ring. The force exerted by MotA 
is modeled as a constant force along an energy surface, 
and the conformational transitions are described as hops 
between the two respective surfaces (Fig. [2]). 

The rotation dynamics of the rotor is given by the fol- 
lowing ovcrdamped Langcvin equation: 

den ^ du^^je,) 

Th^ = -L^^+^l + ^rW (1) 



II. THE STATOR-ROTOR INTERACTION 

The interaction between the stator proteins and the 
rotor proteins is modeled according to the model of Oster 
and coworkers 0, which is based on the description of 
Blair and coworkers [1^1 ■ Since our model builds upon 
this model, we briefly describe its main ingredients. For 
details, we refer to the Supporting Information. 

According to the proposal of Blair and Oster, the mo- 
tor cycle of each stator protein consists of two "half 
strokes" . During the flrst power stroke of a stator pro- 
tein, two protons bind MotB residues of the stator pro- 
tein [l^ (Fig- HI)- This leads to a thermally activated 
conformational transition, which allows one MotA loop 
to exert a force on one FliG protein of the rotor protein 
complex (Fig. [1]). During the second stroke, the recov- 
ery stroke, the two protons are released to the cytoplasm, 
triggering another conformational transition of the sta- 
tor, and allowing another MotA loop to exert a force on 
the FliG protein. At the end of this cycle, the rotor has 



Here, 7r is the friction coefficient of the rotor; Us {Oj) 
is the rotor-stator interaction energy as a function of 
9j = 9n — (Fig. [2]), where 0r denotes the rotor rota- 
tion angle, the angle of the immobile stator protein 
J, and Sj is a binary variable denoting the conforma- 
tional state of stator j; rjn it) is a Gaussian white noise 
term of magnitude Y^SfceTTR; is the number of stator 
proteins. The torque Fl denotes the external load. As 
discussed in [1, [13, EBIj for the system studied here, the 
torque-speed curves under conservative load and viscous 
load are identical. However, as we will show below, the 
type of load does markedly affect the switching dynamics. 

The hopping rate for a stator protein to go from one 
energy surface to another depends upon the free-energy 
barrier separating them. We make the phenomenological 
assumption that the hopping rate depends exponentially 
on the free-energy difference, in a manner that obeys de- 
tailed balance: 

(2) 
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Here, kg sets the basic time scale, and MJs-s'XOj) ~ 
Us'_{Oj) — Usj {Oj). The function w{9j) describes the pro- 
ton hopping windows (Fig. [2]), which reflect the idea 
that the proton channel through the stator is gated by 
the motion of the rotor 0, [3 ■ 

Fig. SI of the Supporting Information shows that this 
model, developed by Oster and coworkers [1], accurately 
describes the torque-speed relation of the flagellar motor 
of E. coli. 



III. THE ROTOR SWITCHING DYNAMICS 

In E. coli, the fraction of time the motor rotates in the 
clockwise direction, the so-called clockwise (CW) bias, is 
controlled by the concentration of the intracellular mes- 
senger protein CheYp. This protein modulates the CW 
bias by binding to the ring of FliM proteins. This ring 
is connected to the ring of FliG proteins, which interact 
with the stator proteins (Fig. [T|). 

The molecular mechanism of the switch is unknown. 
Yet, it is widely believed that the binding of CheYp to 
FliM tends to change the conformation of FliM. Follow- 
ing earlier work, we assume that each FliM protein can 
exist in either a CW or CCW conformational state and 
that binding of CheYp shifts the relative stability of these 
two conformational states [l^, [23|. Moreover, we as- 
sume that also each FliG protein can exist in cither a CW 
or CCW conformational state. In the spirit of an MWC 
model 0, we assume that the energetic cost of having 
two rotor protein molecules in two different conforma- 
tional states is prohibitively large. We can then speak of 
the rotor being in either the CW or the CCW state. 

When the rotor complex switches from one state to 
another, the interactions between the FliG proteins and 
the stator proteins change, due to the new conforma- 
tional state of the FliG proteins. In our model, each 
stator-rotor interaction is described by 4 energy surfaces, 
with the subscript Sj = 0, 1 denoting the conforma- 
tional state of the stator protein j and the superscript 
r = 0, 1 denoting the conformational state of the rotor 
protein (clockwise or counter-clockwise) . We assume that 
the two rotor surfaces corresponding to a given state of 
the stator are simply each other's mirror image (the po- 
tentials are flipped in the 9 direction), but offset by an 
energy difference AC/bias that is given by the CW bias 
Pew = exp(-/3A[/bias)/(l + exp(-/3AC/bias) (Fig. HP). 
This yields the following instantaneous rotor switching 
rate: 

F^^'({0,}) = fcoCxp[AC/'-'-'({%})/2], r,r' = 0,1, (3) 

where AC/-'({0,}) = E^t'i < ) " U^^idj)- Impor- 
tantly, the instantaneous switching rate does not depend 
on the load, although the average, effective switching rate 
does, as discussed below. Note also that our model as- 
sumes that the kinetics of CheYp binding is fast on the 
time scale of switching, in contrast to a recent model Q. 



Fig. [3] shows the switching dynamics when the load 
is conservative. The conservative load is modeled as a 
constant torque in a direction opposite to that of the ro- 
tation direction of the rotor (Fig. [T|); after the rotor has 
switched direction, the conservative force instantaneously 
changes sign. This means that the switching dynamics of 
the motor with a conservative load reflects the switching 
dynamics of the rotor complex; in the next section, we 
study the effect of the flagellum, by studying the switch- 
ing dynamics using a viscous load, which depends on the 
dynamics of the flagellum. 

Fig. [3J^ shows the average switching rate in the for- 
ward and backward switching direction as a function of 
the conservative load. As expected, the CCW CW 
switching rate increases as the CW bias increases. More 
interestingly, the switching rate increases exponentially 
with the external load. As we describe below, this is key 
to understanding the bump in the power spectrum. 

The exponential dependence of the switching rate on 
the load can be understood by noting that the effective 
switching rate is given by 

K^iti^ = / de^p{en)k^-^'i{e,}), (4) 

where P{Or) is the stationary distribution of the ro- 
tor's position. Increasing the load shifts P{0^) to val- 
ues of where the driving force for switching, AU"^^^ , 
is larger (Fig. [2|3). Since the instantaneous switching 
rate k^^^ ({^'j}) depends exponentially on AU^^^ ({^j}) 
(See Eq. [3]), the effective switching rate k^^[^^^ increases 
strongly with load in the low-load regime, as observed in 
the experiments of Fahrner et al. jlOj | . 

These experiments also suggest that in the high torque 
regime, the switching rate decreases with the load jlOj . 
which is not captured by the model presented here; yet, 
we believe that if we would introduce a switching window 
analogous to the proton hopping window introduced by 
Xing et al. to describe the torque-speed relation 0, we 
could reproduce this. However, this is not critical for the 
problem considered here, since the switching dynamics of 
a motor with a flagellum, as described in the next section, 
is mostly determined by how the switching rate changes 
with load in the low to intermediate torque regime, shown 
in Fig. [21 

Fig. [2|3 shows the power spectra of the switching dy- 
namics. It is given by a Lorentzian, which shows that the 
switching of the rotor without a flagellum can be modeled 
as a random telegraph process. 



IV. FLAGELLUM DYNAMICS 

In the model discussed above, after a switching event, 
the torque on the motor immediately changes sign and 
instantaneously reaches its steady-state value. However, 
in the experiments of Cluzel and coworkers, the switching 
of the motor was visualized via a bead that was attached 
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FIG. 3: Switching dynamics with conservative load - the load 
is constant in magnitude, but instantaneously changes sign 
upon a rotation reversal. A: Switching rate as a function of 
the load tl in the forward CW — > CCW (fc+) and backward 
(CCW ^ CW) direction (fc_) for CW bias = 0.1, 0.5, and 0.9. 
Note that due to the symmetry of our model, the switching 
dynamics in the forward (backward) direction for CW bias = 
a;, equals the switching dynamics in the backward (forward) 
direction for CW bias = 1 — x. B: Power spectra S(lo) for CW 
bias = 0.1, 0.5 and 0.9. 



to the flagellar filament [1, 01 ■ We argue that the flagel- 
lum dynamics is critical for understanding the switching 
dynamics of the flagellar motor. 

Darnton and Berg recently studied polymorphic transi- 
tions of a single filament using optical tweezers • The 
following three observations were made: 1) The transi- 
tions occur in discrete, rapid steps that are stochastic in 
nature, suggesting that they are activated processes dur- 
ing which an energy barrier is crossed; 2) In between the 
steps, the filament behaves as a linear elastic object that 
accumulates elastic strain energy that is released during 
the next transformation; 3) During a step, not the whole 
filament is converted, but micrometer-long sections. 

On the basis of these three observations, wc have con- 
structed the filament model shown in Fig. |4l It con- 
sists of a number of harmonic potentials as a function of 
the winding angle 0, corresponding to different confor- 
mational states of the filament. The left-most well cor- 
responds to the normal state, which is the polymorphic 
form of the filament when the motor runs in the CCW 
direction. The right-most well corresponds to the curly 
state, which is one of the polymorphic forms that the 
filament adopts when the motor runs in the CW direc- 
tion. The states in between correspond to an ensemble 
of polymorphic forms that includes not only the coiled 
and semi-coiled states, but, according to observation no. 
3 above, also states in which one section of the filament 
is in one distinct polymorphic state, while the other is 
in an another conformational state. According to obser- 
vation no. 2, and following Darnton and Berg we 
assume that the free energy of a filament in a given 
polymorphic state m is quadratic in the curvature and 
torsion (see Supporting Information). The curvature k 
and torsion t are functions of the height of the filament 
z and the winding angle 9. We assume that at each in- 
stant, the height has relaxed to its steady-state value, 
which means that becomes a quadratic function of 
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FIG. 4: Energy surfaces for the flagellum. The left-most curve 
(m — 1) corresponds to the normal state, the right-most curve 
(m = A'') corresponds to the curly state, while the interme- 
diate conformational states correspond not only to the semi- 
coiled state, but also to hybrid filaments consisting of different 
sections of these polymorphic forrtis, as observed in the pulling 
experiment of Darnton and Berg [14| . The polymorphic tran- 
sitions are modeled as stochastic jumps between these energy 
surfaces (see Eq. [6]). 

the winding angle only: 

Uli9) = lkei0-9r,,f, (5) 

where the torque constant kg is given by the Young's 
and shear moduli and the contour length of the filament; 
the value chosen is consistent with the measurements of 
Block et al [2l| and Darnton and Berg (see Sup- 
porting Information) . For simplicity, we assume that the 
potentials are equally spaced, and have the same torque 
constant and well depth, although under neutral pH the 
normal state is the most stable one [l^l- The total dif- 
ference in winding angle between the normal (left-most) 
and curly (right-most) state is about 50 rounds, which is 
the correct order of magnitude based on the clastic prop- 
erties of the filament (see Supporting Information). This 
is an important parameter, since it directly affects the 
characteristic switching time. 

Motivated by observation 1, we assume that the tran- 
sition from one conformational state to another is an ac- 
tivated process, with a rate constant 

fc„.(0) = keMiUliO) - Ul,{9))/2]. (6) 

Denoting the position of the load (bead) with 6'l, the 
load dynamics is given by 

lL^^-kei9L-9n.-em)+VLit). (7) 

Here, 7l is the friction coefficient of the bead, and t^l is 
a Gaussian white noise term of magnitude y/2kBTj[. 

Fig. [5] shows the switching characteristics of this sys- 
tem. Wc show the dynamics of the bead instead of 
the motor, since that has been measured experimentally; 
however, the switching dynamics of the two are very sim- 
ilar. Fig. [SK shows the distribution of waiting times, 
for CW bias = 0.1, 0.5, and 0.9. These distributions 
agree remarkably well with those observed by Cluzel and 
coworkers Q ■ Firstly, the distributions are not exponen- 
tial, as would be expected for a random telegraph pro- 
cess: The distributions exhibit a clear peak at around 
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FIG. 5: The switching dynamics of a motor with viscous load 
(see Eq. [7}. Panel A: Distribution of waiting times for the 
forward CW — > CCW transition (fc+) and backward CCW 
CW transition (k-), for CW bias = 0.1, 0.5 and 0.9. Panel B: 
The power spectra S{uj) for CW bias = 0.1, 0.5 and 0.9. Our 
model is symmetric by construction — the CW energy surface 
is the mirror image of the CCW surface (Fig. [2)3) and the 
wells of the filament potential are of equal depth (Fig. |4|. 
Accordingly, the distribution of the forward (backward) tran- 
sition for CW bias = x overlaps with that of the backward 
(forward) transition for CW bias — 1 — x. 



0.4 s~^. Secondly, the waiting-time distribution for the 
forward (CW-^CCW) transition shifts from a narrow dis- 
tribution at CW bias = 0.1 to a broad distribution at CW 
bias = 0.9. Moreover, the position of the maximum of 
the distribution shifts to longer times. All these features 
are in near quantitative agreement with experiment. 

Fig. [5j3 shows the power spectra of our model, for CW 
bias = 0.5, and for CW bias = 0.1, 0.9 (they are identical 
because of the symmetry of our model) . It is seen that the 
spectra exhibit a distinct peak at w ^ 1 s~^. Moreover, 
the peak is most pronounced when the CW bias = 0.5. 
These characteristics are observed experimentally Q ■ 

Our model predicts that the peak in the power spec- 
trum arises from the interplay between the conforma- 
tional dynamics of the flagellum and the dependence of 
the switching rate on the load (Fig. [Sj'V). The idea is illus- 
trated in Fig. [51 After a switching event of the rotor, the 
torque is initially in the original direction, but decreases 
rapidly in magnitude as the filament reaches its optimum 
winding angle (Fig. [U Fig. [BJ\); in this regime, the load 
on the motor is negative, and the elastic strain energy in 
the filament decreases. As the rotor drives the filament 
beyond its optimal winding angle, the torque changes di- 
rection and increases in magnitude; the load on the motor 
becomes positive, and the strain energy in the filament 
builds up. This strain energy can, however, be released 
via a polymorphic transition, leading to a sudden change 
in the direction of the torque. This process repeats it- 
self until the filament reaches its final polymorphic form, 
upon which the strain energy can no longer be released, 
and the torque increases to reach a plateau when the 
viscous drag equals the motor torque. The peak in the 
power spectrum can now be understood by combining the 
time trace of the torque (Fig. [^^^-) with the dependence 
of the switching rate on the load (Fig. [3]A), yielding the 
switching propensity (i.e., the probability to switch per 
unit amount of time at a given time t) as a function of 
time, as shown in Fig. [BJ3. After a switching event, the 



FIG. 6: The mechanism of switching. A: Typical time traces 
for the torque, motor and bead position, for CW bias — 0.5. 
B: The switching propensity as a function of time after a 
switching event, for CW bias = 0.1, 0.5, and 0.9. This figure is 
obtained by combining the switching propensity as a function 
of the conservative load, as shown in Fig. [3}\., with the average 
force as a function of time after a switching event, as shown 
in panel A. 



torque flip-flops around zero and the switching propen- 
sity remains low. But when the flagellum has reached 
its final polymorphic form, the torque can no longer be 
released, and the switching propensity increases signif- 
icantly. The peak in the spectrum is precisely caused 
by the fact that the switching propensity function is not 
constant in time, as for a Markovian Poisson process, but 
increases with time. 



V. COARSE-GRAINED MODEL 

Our calculations suggest that a useful coarse-grained 
model for understanding the switching dynamics is one in 
which the system stochastically flips between two states 
with time-dependent propensity fimctions (Fig. [6j3): 



fe+(t) 
CW ;=i CCW, 

fc_(i) 



(8) 



where the propensity functions are given by the following 
piece- wise linear functions: 



+ (fcr'^ - kfn^-^ h<t<t2 (10) 

t2 — ti 



t < ti (9) 
10) 

t > t2 (11) 



The important parameters of this model are the lag time, 
Ta = (ti-l-i2)/2, the minimum and maximum propensity, 
/c™'" and fc™'*'^, respectively, and to a lesser extent the 
sharpness of the transition = [k'^^'^ — A:™'")/(t2 — ^i)- 
For this model, the waiting-time distribution and power 
spectrum can be obtained analytically (see Supporting 
Information) . 

The maximum propensity function k™^^ depends on 
the CW bias, which determines how the rotor switch- 
ing propensity depends upon the load (Fig. [3lA), and 
the maximum load itself, which is determined by the 
drag coefficient of the load, 7l, and the torque-speed re- 
lation; to a good approximation, the maximum torque 
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is given by the intersection of 7l times the speed (the 
load hne) and the torque-speed curve The mini- 

mum propensity function fc™™ does not depend upon this 
maximum load, but rather upon the torque at which the 
flagellum undergoes a polymorphic transition — the poly- 
morphic transitions release the elastic strain energy be- 
fore this maximum load is reached, keeping the load and 
hence fc™'" low. The peak in the power spectrum arises 
when fc™'" < k^'^^, which leads to the prediction that 
the peak may disappear when the number of stators is 
reduced. The position of the peak is determined by T^, 
which is given by the difference in winding angle between 
the normal and curly state divided by the average speed 
at which the rotor drives the systems between these two 
states; this speed depends upon the torque-speed rela- 
tion and the load as a function of time. Interestingly, 
polymorphic transitions of filaments of swimming bacte- 
ria occur on time scales of 0.1 s close to the peak of 
the waiting-time distribution, supporting our idea that 
they set the characteristic switching time. 

The dependence of fc™^'^ on the CW bias can explain 
the change in the waiting-time distributions when the 
CW bias is varied. (Fig. [SJ^). When the CW bias is large, 
^max ig^j-gQ^ because the switching propensity at the 
maximum load as set by the balance of the drag and the 
motor torque is large (Fig. [SK)- Consequently, the rotor 
typically switches to the CW state before the switching 
propensity can reach its plateau value. This explains the 
narrow distribution of CCW intervals when the CW bias 
is large, as observed in both the model (Fig. and 
experiment 0] . For the reverse transition the situation is 
qualitatively different. When the CW bias is large, k'i^^^ 
is low, which means that the system can enter the regime 
in which the switching propensity is constant before it 
switches to the CCW state. This constant propensity 
leads to an exponential tail in the distribution of CW 
(CCW) intervals when the CW (CCW) bias is large, as 
observed in both the distributions of the model (Fig. [SIA) 
and those measured experimentally 



VI. DISCUSSION 

We have presented a statistical-mechanical model that 
describes the switching dynamics of a rotary flagellar mo- 
tor. Its foundation is the assumption that the rotor pro- 
tein complex can exist in two conformational states cor- 
responding to the two respective rotation directions, and 
that switching between these states depends on interac- 
tions with the stator proteins, which also drive the ro- 
tation of the rotor complex. This naturally couples the 
switching dynamics to the rotation dynamics. The load 
docs not directly change the relative stability of the ro- 
tor's conformational states, but it does change how often 
the stator proteins during their motor cycle favor one 
conformational state of the rotor over the other. This, 
according to our model, is the principal mechanism that 
makes the switch sensitive to torque and speed. Another 
key element of our model is that after a switch, it takes 
time for the load to build up, due to the polymorphic 
transitions of the filament. 

Several predictions emerge from our model that could 
be tested experimentally. One is that the change in the 
torque on the filament upon a motor reversal leads to a 
series of polymorphic transitions, which could be tested 
by applying a torque on a single filament using magnetic 
tweezers. Another is that the switching dynamics of the 
rotor without the viscous load of the flagellum is that 
of a two-state Poisson process, in contrast to a recent 
model [1]; moreover, in the zero-load regime, the switch- 
ing propensity increases with the load. These predictions 
could be tested by measuring the rotation dynamics of a 
bead that is connected either directly to the stub, or to 
a very short filament. Our model predicts that its power 
spectrum does not have a peak. 
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